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Abstract. Loss is a critical roadblock to achieving photonic quantum-enhanced 
technologies. We explore a modular platform for implementing integrated 
photonics experiments and consider the effects of loss at different stages of 
these experiments, including state preparation, manipulation and measurement. 
We frame our discussion mainly in the context of quantum sensing and focus 
particularly on the use of loss-tolerant Holland-Burnett states for optical phase 
estimation. In particular, we discuss spontaneous four-wave mixing in standard 
birefringent fibre as a source of pure, heralded single photons and present 
methods of optimising such sources. We also outline a route to programmable 
circuits which allow the control of photonic interactions even in the presence of 
fabrication imperfections and describe a ratiometric characterisation method for 
beam splitters which allows the characterisation of complex circuits without the 
need for full process tomography. Finally, we present a framework for performing 
state tomography on heralded states using lossy measurement devices. This is 
motivated by a calculation of the effects of fabrication imperfections on precision 
measurement using Holland-Burnett states. 
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1. Introduction 

The field of integrated pliotonics is a promising area for the development of quantum- 
enhanced technologies with applications being pursued in communications, metrology, 
simulation and information processing [1] [21 121 HI 13 E] ■ These quantum experiments 
can be broadly broken up into three stages: state preparation, manipulation, and 
measurement. Over the last decade or so, integrated approaches have made substantial 
advances in each of these three areas, particularly with regard to fibre- and waveguide- 
based nonclassical fight sources [71iaini[Tni[ni[Il[Il[IllI3IIllI71[IHl[Ill2nilII], 
on-chip quantum circuits [lllll[Sl[Tl[21[l[lSl|ll|17],and integrated detection ^51 

Perhaps the greatest challenge facing photonic quantum applications is the issue of 
photon loss, which is generally present in all three of these stages. Loss is particularly 
critical in quantum information processing (QIP) and quantum-enhanced metrology, 
where strict thresholds exist for the efficiencies required for a device to out-perform 
its classical counterpart |3TJ [321 [331 IM] • The most significant contributions to loss 
in photonic experiments generally arise at the interfaces between different stages of 
the experiment, e.g., between sources, which may be either spatially or spectrally 
multimode, and circuits, which often rely critically on low-distinguishability quantum 
interference. In this example, bulk optics experiments typically require the use of 
strong, inherently lossy spatial and spectral filtering to reach the necessary operation 
fidelities. 

The key advantage of integrated optics in overcoming these losses is that it 
provides the potential for complete spatial and spectral control of the underlying 
optical field modes. The tight confinement and guided nature of these modes 
provides exquisite control of their spatial properties, both in terms of their transverse 
profiles and paths of propagation. This enables strong interactions between highly 
indistinguishable optical modes. It is therefore possible to fabricate devices with 
precise, stable, and potentially complex circuit configurations which, because they are 
small and monolithic, are inherently stable without the need for complex stabilisation 
techniques and provide high- visibility interference [35l [U [211 [21 [36]. This holds 
great promise for many photonic QIP applications, which can require complex 
circuits containing many nested interferometers [37l |38l [39l |40l [41] , by alleviating the 
impractical space and stability requirements facing bulk-optics implementations. At 
the same time, by carefully managing the dispersion properties of integrated devices, 
either by natural means (material selection) or engineering design (e.g., poling or mode 
confinement), it is also possible to exercise fine control over the spectral properties of 
the optical modes. Such control can be used, for example, to design spectrally pure, 
heralded single-photon sources [H [42] [43l [20l [21] , a critical requirement for scaleable 
photonic QIP systems. 

The holy grail of integrated photonics is a single device containing an entire 
photonic quantum experiment, state preparation, manipulation, and measurement. 
Such an approach, however, faces an enormous challenge; a fully integrated platform 
requires the marriage of many fundamentally different, often-incompatible technologies 
within a single material substrate. For example, thermo-optic phase control [H [2] 
would clearly be impossible on a device which is being cooled to cryogenic temperatures 
(mK) in order to operate a superconducting TES (transition-edge sensor), which 
currently provide the leading performance in terms of both quantum efficiency and 
photon-number resolution. 
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An alternative approach, achievable within the current state of the art, is to 
employ a more modular solution, where each stage of the experiment is performed 
using the technology platform that is most suited to the task. The main task of 
integration is then simply to ensure that the interface between platforms does not 
introduce an unacceptable level of loss. This is not an unreasonable target, however, 
since the unprecedented spatial-mode control available in integrated architectures 
already provides the means for optimising the coupling efficiencies between platforms. 
Such an approach also presents a strong advantage in terms of flexibility. This cluster 
of technologies would become a platform on which large scale complex experiments 
could be performed, enabling progress in both practical applications of quantum 
enhanced technologies and fundamental research, and could include components 
for active phase control [21 E], Bragg reflectors for frequency filtering and building 
cavities |42| , and on-chip micro-fluidic capillaries for interfacing photonic systems with 
fluids or fluid suspensions [H], to name just a few examples. 

In this paper we will focus on the specific example of sensing in a real-world 
integrated photonic device via (linear) quantum interferometry as it provides a clear 
example of the stringent requirements placed upon the different modules. We will 
follow the chain of state generation, manipulation, and detection, discussing in turn 
the key requirements for and issues surrounding each. We will also introduce a 
technique for characterising optical elements in an integrated photonic device, as well 
as a framework which allows tomographic reconstruction of a heralded state, including 
all photon number subspaces. 



2. Quantum sensing — an example context 
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Figure 1. Schematic of a typical metrology scenario. (a) General two 
mode interferometer with inefficient state preparation, two lossy modes, and 
inefficient detection, (b) A detailed equivalent network involving a Mach-Zehnder 
Interferometer (MZI) with the losses specified by individual beam splitters of 
transmissivity rjp, r), and rj^. 



Quantum-enhanced metrology is concerned with a single task: preparing a 
quantum state which is sensitive to a parameter, ^, and implementing a measurement 
on that state so that the uncertainty in the measurement of the parameter, A(/), is 
lower than the uncertainty that would be obtained using the same number of classical 
resources. Typically, only the scaling behaviour of with the number of particles in 
the state, iV, is considered since for large enough N a class of states which provides 
only a constant factor improvement can always be beaten by one which has a more 
favourable scaling. The tightest known limit on precision is given by the Heisenberg 
limit (HL), A0 > /S.(f>HL = ^/y/vN , where v is the number of experimental trials. 

The situation we will consider in this paper is shown schematically in figure [Ija). 
For any given state and measurement, the input state p evolves to become p(4>) 
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and a measurement is performed with outcomes represented by the operators {H^}. 
Performing v trials provides a data set and the precision with which can be 
estimated from this data set, Acj), is bounded by the Cramer-Rao bound (CRB), 
A(/) > l\(j)cRB — ^1 \fvF{^- Here, is the Fisher information (FI), a measure 

of the information which can be obtained about from the data set. The CRB 
can be saturated for large v by maximum likelihood estimation. Interestingly, a 
measurement-independent version of the FI can be found by maximising over all 
physical measurements giving the quantum Fisher Information (QFI), Fqi^c^). This, 
in turn leads to the lowest possible bound on the precision obtainable using a given 
state, the quantum Cramer-Rao bound (QCRB) A0 > Ac^qcrb = 45 . 

To judge whether or not a device could perform better than its classical 
counterpart, we must compare the QCRB of the prepared state to that of N completely 
uncorrelated photons. In the ideal situation of a perfectly transmissive device, this 
leads to the standard quantum limit A(/) > A(f>sQL = \l\JvN . In the presence of loss, 
however, this becomes the standard interferometric limit (SIL) [35], > A(j)siL = 
1 / y/urjrjdN where 77 and 77^ model transmission inside the interferometer, and detection 
efficiency respectively (see figure [T]) . The SIL is unaffected by the parameter rjp as 
a state which saturates the SIL, the phase-averaged coherent state, is unaffected 
by rjp apart from a scaling in average photon number. It is therefore considered 
that a phase-averaged coherent state with arbitrary average photon number can be 
prepared perfectly. In general, however, all losses including rjp play a crucial role in 
the QCRB for any given prepared state. If Ac^qcrb < AcpsiL then the state which is 
prepared could perform better than its classical counterpart, but it will only do so if 
an appropriate measurement scheme is implemented. If such a scheme is implemented, 
A(pcRB < Acj)siL and the device truly does outperform its classical counterpart. 

3. State generation 

The most widely studied state for quantum-enhanced metrology is the NOON 
state [421 EHl SSI m El EHl El]. This path-entangled state is a superposition of N (0) 
photons in one mode with (N) photons in the other and, under perfect transmission, 
it saturates the HL. It does, however, suffer from a major drawback in the lossy 
case. When even a single photon is lost from the NOON state, the remaining state is 
entirely insensitive to the phase 0. Thus all photons must make it through the device 
and be detected in order for information to be gained about </>, an exponentially 
unlikely event for increasing N [33]. Schemes for generating high-A^ NOON states 
are also complex, often requiring many nonlinear elements [52] or carefully aligned 
cavities [53] . Optimally loss-tolerant A^-photon states have also been studied [46l [54] , 
but preparation strategies and detection schemes for these states have not yet been 
found and knowledge of the exact channel loss would also be required. 

By contrast to NOON states and optimal states, Holland-Burnett (HB) or twin- 
Fock states [5S] provide a much simpler route to quantum-enhanced metrology whilst 
retaining near-optimal loss tolerance. An HB(A^) state is prepared by interfering two 
A^-photon Fock states at a 50:50 beam splitter (see figure [TJ)). The resulting state 
retains phase sensitivity under loss both before and after this beam splitter and hence 
does not need perfect transmission to be useful. As well as this, the QFI for the HB(A^) 
state is attained through the use of photon-number-resolving detectors, regardless of 
the loss parameters [34]. The requisite Fock states can be prepared by heralding 
A'^ photons from two two-mode photon-number-correlated sources, for example based 
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on spontaneous parametric downconversion (SPDC) or spontaneous four-wave mixing 
(SFWM). Recent work [HJ has put bounds on the rjp, rj, and 77^; that are necessary to 
beat the classical limit, showing that for 77^; = 0.6 and rj = 0.95, a state preparation 
efficiency of rjp > 0.91 is required. Even for the best currently available detectors with 
rjd = 0.98, the requirement on rjp is relaxed only to rjp > 0.71. This puts extremely 
challenging bounds on the current state of the art for quantum sources. 

The key ingredient for generating HB(A^) states is the creation of A^-photon Fock 
states with sufficient purity to achieve the quality of quantum interference required 
for high-fidelity state production. More generally, if these sources are to be useful 
in a scaleable way for integrated photonics quantum technologies, they will need to 
be either on-demand or heralded. A heralded source can be converted into an on- 
demand source using a controllable quantum memory with a high time-bandwidth 
product [SSI [SH] . The important characteristics of such sources are heralding or 
preparation efficiency, purity, and indistinguishability of the photons produced (in 
polarization, spatial and temporal degrees of freedom) . In addition, in order to avoid 
excess loss at the interface with the quantum circuit, the photon mode structure 
should be compatible with integrated devices, particularly in terms of spatial overlap, 
to ensure optimal performance. If the photons produced by these sources do not 
adhere to strict requirements for these parameters, the performance of the quantum 
protocols rapidly decays (e.g., [50] ). 

One of the most common and most successful ways to produce heralded photonic 
quantum states makes use of nonlinear processes such as either spontaneous parametric 
down-conversion (SPDC) or spontaneous four- wave mixing (SFWM). These are 
robust, room-temperature processes that can produce large numbers of photons for 
moderate pump powers. Because of the strong photon-number correlations induced 
by the nonlinear interaction, they can be used to produce high-quality heralded Fock 
states. Unfortunately, however, the generated state also contains terms with more 
than the desired number of photon pairs in a single pulse. This ambiguity leads to 
dramatically reduced performance in many QIP applications if standard avalanche 
photodiodes (APDs) are used as heralding detectors, because they are unable to 
distinguish between different photon- number terms |60| . In the ideal case, this 
problem can be largely mitigated by instead using photon-number-resolving heralding 
detectors [29l[61]. In the presence of loss, however, it becomes exacerbated and cannot 
simply be overcome by using a stronger pump to increase the count rate. Loss acts to 
weaken the photon-number correlations in the generated state and therefore decrease 
the quality of interference exhibited by the heralded Fock states. Furthermore, 
although increasing the pump does increase the signal strength, it can increase the 
strength of the noise-contributing higher-order terms by a larger proportion, because 
although they occur less frequently than the signal term, they are often more likely to 
be detected when they do occur. Thus, quantum state production must be optimised 
to provide the highest possible heralding efficiency. 

In general, photon pair sources based on SPDC or SFWM produce pairs which 
display spectral entanglement and this manifests itself in a joint spectral amplitude 
(JSA) which is not factorable. A typical example of such a JSA is shown in figure [2ja) 
where the correlations are clearly visible as correlations in the frequencies of the 
two photons. Consequently, although detecting N photons in one mode heralds the 
presence of N photons in the other, the heralding detectors are typically spectrally 
non-selective and therefore project these heralded photons into a spectrally mixed 
state. The usual technique for removing these spectral correlations is to use narrow 
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spectral filters, but this improvement in purity comes at the expense of decreased 
heralding efRciency, i.e. decreased rjp, which therefore also increases the effects of 
noise in the experiment. 

The trade-off between heralding efficiency and spectral purity for the highly 
correlated joint spectral amplitude of figure [Sf a) is shown in figure [2ljb). The purity 
of the unfiltered heralded signal photon is 27%. Using spectral filters that only pass a 
symmetric portion of the joint spectrum dramatically decreases the heralding efficiency 
of the source so that, for a purity of 95%, heralding efficiency is limited to 9.7%. 
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Figure 2. (a) A typically correlated joint spectrum for a spontaneous four- 
wave mixing (SFWM) photon-pair source, (b) Heralding efficiency (solid) and 
purity (dashed) as a function of Gaussian filter bandwidth (in units of the pump 
bandwidth) for the correlated case, (c) Joint spectrum for a SFWM photon- 
pair source with matched phase matching and pump bandwidths. (d) Heralding 
efficiency (solid) and purity (dashed) as a function of Gaussian filter bandwidth 
(in units of the pump bandwidth) for the matched bandwidths case. 



A promising route to generating heralded pure-state single photons with high 
efficiency and purity is to directly engineering the nonlinear optical interaction [621 
I63| I64| [13, 65, 66^ 167) . In these approaches, an intense pump is converted into a pair 
of photons, denoted signal (s) and idler (i) mediated by the nonlinear medium. This 
process must obey both energy and momentum conservation, which generally leads 
to spectral and momentum correlations. However, by carefully choosing material 
parameters and pump configuration, these correlations can be eliminated. SFWM 
utilizing birefringent phase matching in optical waveguides offers several advantages 
over other nonlinear optical sources, such as perfect spatial mode matching with 
existing photonic circuitry allowing low-loss integration, flexible tuning of the signal 
and idler central wavelengths through birefringence, and tolerance to waveguide 
properties tl8j,68j. In this scenario, the pump propagates along the fast axis while the 
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signal and idler photons are produced along the orthogonal slow axis of the waveguide. 
The wave- vector mismatch includes a term proportional to the birefringence and pump 
central frequency. Thus, by adjusting the pump central frequency and birefringence 
one can tune the central frequencies of the signal and idler photons. To achieve 
the nearly factorable states required for heralding pure-state single photons with 
birefringent phase matching, the pump bandwidth must be balanced with the phase 
matching bandwidth which is determined by the interaction length |181 168j . A typical 
JSA produced under such balanced conditions is shown in figure [2jc). The correlations 
are significantly reduced despite the remaining lobes which result from the sinc-like 
form of the phase-matching function. Figure [2jd) shows the effect on the heralding 
efficiency and purity of applying spectral filters to this JSA. For a purity of 95%, 
heralding efficiency in this case is greater than 80% directly demonstrating the higher 
intrinsic purity of the source. 

4. Manipulation 

Once pure single photons have been created, they must be manipulated in such a way 
as to implement, for example, logic gates or a metrology scheme. Only two components 
are necessary in order to implement all unitary operations: beam splitters and phase 
shifters. Beam splitters can be constructed using evanescent [Ml IS], x-couplers [T], 
or multimode interference devices |26| while phase shifters can be implemented by 
thermo-optic [U [2] or electro-optic [69] means. 

Evanescent couplers are fabricated by bringing together two waveguides so that 
the evanescent fields of the two modes overlap, creating a coupling between them. The 
field then flips back and forth between the modes as they propagate in parallel. The 
rate of this flipping is determined by the overlap of the two modes which is dependent 
on the separation of the waveguides. For such a device with a given waveguide 
separation, the splitting ratio is determined by the length of the coupling region, and 
the wavelength being used. They can be produced using either lithographic or direct- 
write techniques. X-couplers are produced by physically crossing the waveguides at a 
shallow angle giving a splitting ratio which is a function of this angle. Direct write 
techniques must be used to fabricate x-couplers since in order to produce this shallow 
angle, sharp features in the guiding structure are required which cannot currently be 
produced using lithographic techniques. 

A detailed comparison of thermo-optic phase shifters against electro-optic versions 
is yet to be performed for quantum devices because, although electro-optic phase 
shifters are common in commercially available classical integrated devices, they have 
not been demonstrated in an integrated photonic device. They have several potential 
advantages over their thermo-optic counterparts, however, most notably a reduced 
switching time which, if combined with integrated quantum memories, could provide 
a route to one-way photonic QIP. 

For the specific example considered in this paper, it is critical to be able to produce 
beam splitters with accurate splitting ratios. The splitting ratios of the two beam 
splitters shown in figure[Hb) affect both the fidelity of the HB state produced and the 
measurement scheme which is implemented, limiting the FI that can be obtained. 

In order to show this, we begin by defining the angular momentum operators in 
the Schwinger representation as: 

J. = ^(a^fo + aSt), = l(at6-aSt), J, ^ ^{a^a -b^b), (1) 
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J2^J2^J2 + J2. (2) 

In a basis defined by the mutual eigenbasis of and J^, the ideal input state \N)a\N)b 
appears as \N, 0). In this picture, the phase operator is |7D] 

P{<l>)=e^^'% (3) 

while a beamsplitter is given by 

= 6-2'"^% (4) 

where sin^ 6 = v is the input beamsplitter transmissivity. The effective phase operator 
is then given by 

P{(j)) = B{e)P{(t))B{9)^ ^ g.0(cos2ej,-si„2ej„)^ 

where we have used the SU(2) commutation relations of the angular momentum 
operators. For a 50/50 beamsplitter, 9 = 7r/4, in which case the effective role of 
the beamsplitter and the phase combined is to effect a rotation about the —Jy axis. 
In general, this axis of rotation is given by ([S]), and given that the input state points 
along the z-axis, the radius of rotation shrinks by a factor of sin 

as depicted in figure |3Ja). The net effect is to suppress the phase picked up by the 
state \N)a\N)i, with the beamsplitter ratio v. Since the phase picked up is suppressed 
by this factor, the precision is also suppressed by the same amount so that the QFI 
in this case is given by 

Fq^4:N{N + l)v{l-v). (6) 

This shows that the best attainable precision is provided by a 50/50 input 
beamsplitter. In that case, we know that this can be attained with another 50/50 
beamsplitter at the other end of the interferometer and an |A^)(A^| (g) |iV)(iV| detection 
(in the loss-less case). For every other situation, however, there remains the issue of 
attaining the reduced bound provided by ([B]). Allowing for a general beamsplitter at 
the detection end makes the calculation somewhat involved, so we restrict ourselves 
to the |l)a|l)b case. Then the classical FI, F, as a function of the input beamsplitter 
V and output beamsplitter w is shown in figure |3l[b) , assuming that we have perfect 
photon-number resolving detectors. 

In practice, fabrication tolerances result in beamsplitters which do not have the 
reflectivity they were designed to. Fortunately, the MZI provides a solution, since the 
combination of two beam splitters and a phase shifter can be thought of as a single 
beam splitter whose reflectivity, Wcff is a function of the applied phase, 

VcS = vw -f (1 — v){l ~ w) + 2^v{l ~ v)w{l — w) cos (j). (7) 

If both beam splitters were to have a splitting ratio of 1/2 then, by adjusting the 
phase, the full range of reflectivities from perfectly reflective to perfectly transmissive 
can be realised. As the reflectivities of the beam splitters vary away from 1/2, the 
effective reflectivities which can be reached by tuning the phase is reduced, however 
this technique dramatically increases the robustness of the system to fabrication 
tolerances. Figure S] shows the values of the two beam splitters for which reflectivities 
of 1/2 and 1/3 can be obtained. In principle, one could even combine four nominally 
50:50 beam splitters with three phase shifters to produce a single programmable beam 
splitter which is tunable over the full range < Woff < 1. 

Since variation in beam splitter reflectivities is inevitable, a way of accurately 
characterising them is necessary in order to calibrate and predict the behaviour of 
any device. Full process tomography becomes extremely complex and unwieldy for 
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Figure 3. (a) Depiction of the effect of varying v and w on the effective phase 
operator, (b) The Fisher information of the state generated by |l)a|l)i, impinging 
on a beam splitter of reflectivity v when a detection event |l)(f | 5? |1)(1| occurs 
after a beam splitter of reflectivity w, maximised over the phase <j>- 
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Figure 4. (a) The values of v,w for which an effective reflectivity of 1/2 can 
be reached, (b) The values of v,w for which an effective reflectivity of 1/3 can 
be reached. In both plots the reflectivity combinations in the white areas do not 
allow tuning to the target effective reflectivity. The phase required to reach the 
target effective reflectivity is represented by the colour. 



large systems so it is necessary to characterise components individually to build up the 
characteristics of the device. Here we present a ratiometric characterisation technique 
which allows measurement of the splitting ratio independently of the input and output 
coupling efficiencies associated with launching light into and collecting light from the 
device. The situation being considered is shown schematically in figure [51 Although 
only a single beam splitter is shown, this technique can be applied to any beam splitter 
embedded in a larger circuit as long as both inputs to and both outputs from the beam 
splitter can be accessed independently. 

Initially, coherent light with intensity / is coupled into mode a. This light 
propagates through the circuit and is detected at the outputs, having experienced 
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Figure 5. Schematic of a device consisting of a single beam splitter of 
reflectivity u with input and output coupling efficiencies modelled by beam 
splitters with ria,rii, for input coupling transmissivity, and r^cVd ^or output 
coupling transmissivity. The spatial modes are labelled a, b, c, d. Detection 
is performed by photodiodes (PD). 



propagation, coupling, and detection losses. These can be represented by a single 
transmissivity at the inputs, rja and rji,, and a single transmissivity at the outputs, 
77c and rjd.. At the outputs, two intensities are measured, lac and lad in modes c and 
d respectively. These intensities can be written as a function of the transmissivities, 
beam-splitter reflectivity, and the initial intensity, 

lac = VaUVcI, (8) 
lad = Vai^ - u)l]dl ■ (9) 

Similarly, the same light is coupled into mode b, giving 

he = VbO- - u)r]cl (10) 

hd = riburidl. (11) 

Finally, by taking the ratio r = ladhc/ lachd and solving for u, taking the solution 
which gives < u < 1, we get 

u^—^. (12) 




0.05 0.1 0.15 0.2 
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Figure 6. The cff'cctivc reflectivity of a MZI as a function of the power dissipated 
by a thermo-electric phase shifter. Characterisation was performed using a 
ratiometric technique which is independent of coupling losses (see text). 
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Figure [S] shows the resuhs of this technique being apphed to a MZI which is acting 
as a single, programmable beam splitter. The maximum beam splitter reflectivity 
reached is 0.993. The range of applied powers is limited in this case due to the 
robustness of the particular device used. Fitting equation (O to a full fringe would 
allow the individual beam splitter reflectivities to be obtained, however, this equation 
is symmetric under exchange of the two reflectivities and so, although the values can 
be obtained, the ordering is ambiguous. In order to resolve this ambiguity, we must be 
able to measure directly one of the beam splitter ratios. As well as this, the individual 
inputs and outputs of some beam splitters in QIP circuits are not independently 
accessible, for example in a CNOT gate where only three beam splitters may be 
characterised in this way [36] . 

This problem is alleviated by using a camera to collect light scattered out of the 
mode in the transverse direction. The outputs of a particular beam splitter can then 
always be independently accessed, drastically increasing the applicability of the above 
technique. In this situation, the light observed at the camera is a proportion of the light 
propagating in the mode. As with the ratiometric technique outlined above for output 
intensities measured with a linear detector, the proportion scattered does not need to 
be constant for different paths since the ratiometric analysis gives a reflectivity which 
is independent of the scattering ratio and coupling efflciency. Consequently, however, 
the camera used must display a linear response and the image of the output modes 
must not change position on the camera sensor so that the pixel efflciency is constant 
when coupling into each input. 

5. Measurement and Detection 

The final stage in any photonic experiment is detection. High-quality, quantum-limited 
measurements are an essential requirement for demonstrating explicitly non-classical 
behaviour in a quantum experiment, and this is characterised by the detectors, the 
measurement apparatus as a whole, and by the way that apparatus is utilised in the 
experiment. Quantum-limited detection in photonics experiments is normally carried 
out using photon counting. As discussed above in relation to state preparation, in 
the metrology example, the detection must be very efficient if the quantum device is 
to out-perform its classical counterpart in terms of overall precision [331 [M| ■ Indeed, 
even with detector efficiencies close to 100%, the constraints placed on preparation 
efficiencies and propagation losses for achieving true quantum-enhanced operation are 
extremely stringent and have not been demonstrated by any existing experiment. 

The standard workhorse of photon-counting experiments is the APD, a binary 
("zero/many") detector which is readily commercially available. Unfortunately, 
typical quantum efficiencies for silicon APDs reach at maximum up to 60% 
for wavelengths around 800 nm, where they have a peak in sensitivity. By 
contrast, efficiencies of up to 98% [7T] can be reached with the recently developed 
superconducting transition-edge sensors (TES detectors), with the additional 
advantage that they also provide direct photon-number resolution. Both APDs and 
TES detectors are generally operated in a fibre-coupled configuration, making them 
well-suited to a modular integrated-photonics platform. A major drawback with 
the still-emerging TES technology is that the detectors currently have to operate at 
ultra-cold cryogenic temperatures (typically mK). While the efficiency is still limited, 
another way to achieve many of the benefits of photon-number resolution is via 
time multiplexing [29l [61] . Time- multiplexed detectors (TMDs) are pseudo-photon- 
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number-resolving detectors which utihse standard, commercial components and can 
be operated straightforwardly at room temperature. They work by dividing the input 
signal into a number of distinguishable temporal modes, which are then monitored by 
standard APDs. The scaleability of this technique is largely limited by the dead-time 
of the APDs, which is typically around 50 ns for current APDs. This currently leads 
to large readout-response times and therefore limits the repetition rate at which they 
are able to operate. TMDs only provide pseudo-photon-number resolution, because 
there is still some probability that more than one photon will arrive at an APD in a 
single temporal mode — they only provide unambiguous photon-number resolution in 
the limit of many more temporal modes than incoming photons. These effects can be 
in many cases mitigated, however, by appropriately characterising the behaviour of 
the detectors in advance via detector tomography [751 [721 IZl] • 

Detector tomography is a general technique which uses weak-coherent input states 
to probe the detector operation and reconstruct the individual measurement operators 
representing each possible detector outcome, including their description in the photon- 
number basis. To date, this technique has been used to characterise the behaviour 
of both a standard APD and a TMD [72l [73]. An important aspect of detector 
tomography is that it provides a direct estimate of detection efficiency [74] which, as 
we have already discussed, is key to understanding and interpreting observations of 
quantum effects in experiments (e.g. metrology), particularly in the presence of loss. 

A common task that is often necessary in the building and characterisation of 
photonic circuits for quantum applications is quantum state tomography [75l [76l [77] . 
In order to understand the performance of the quantum device, it is important to be 
able to estimate the form of the input quantum state upon which that performance 
depends. This is equally true for integrated photonics applications. Traditionally, 
photonic experiments rely on post-selection, where the state is assumed to be within 
a specific photon-number subspace. This means that loss in the measurement device 
used to perform state tomography can, for the most part, be ignored. Unfortunately, 
however, this is not sufficient in many situations, such as the metrology example 
that we have considered here. In this case, phase information can generally be 
acquired in some degree from all photon-number subspaces, particularly in the case of 
"classical" coherent input states, against which any quantum device must be measured 
to determine whether it improves upon the classical limit. The ability to acquire some 
level of phase information from lower photon-number terms is a signature of loss- 
tolerant probe states, such as the optimally loss-tolerant state [46l [54] and the HB 
states [5S]. Unusually, this is not the case with the NOON state, which is an indication 
of its fragility to loss. However, because of the effects of loss, even in this case it is 
critical to characterise the complete input state, including the contributions from all 
photon-number subspaces, in order to obtain a genuine estimate of the performance of 
the quantum device. Once the complete density matrix of the input state is known, this 
can be used to calculate the QCRB, which defines the usefulness of the prepared state 
for phase estimation. In order to reconstruct the complete density matrix of a quantum 
state, two key experimental ingredients are required; firstly, a heralded source of 
input states to probe the contribution from the vacuum subspace, and secondly, a 
measurement apparatus which provides some level of photon-number resolution, so 
that different Fock subspaces can be interrogated. 

In general, a state tomography experiment can be described by a set of 
detection apparatus settings, labelled here with a, each of which has associated 
with it a set of outcomes, labelled here with 7 so that NaN^ measurements 
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Figure 7. A typical setup used to perform state tomography on p. p is assumed 
to be a single polarisation photonic state in two spatial modes. A MZI and two 
phase shifters allow a tomograpically complete set of measurements to be made. 
Detection is performed by avalanche photodiodes (D1-D3). The spatial modes 
which lead to detectors are labelled mi-m^ while loss modes are labelled lu-lfn 
and ho-ho- 



are implemented. In photonic experiments, the outcomes are usually a particular 
combination of physical detectors firing simultaneously ("in coincidence"). The 
measurements implemented are modelled by a set of positive operator value measure 
(POVM) elements In this scheme, the outcomes for each setting are complete 

so that for every a, J2-y^i — 1- The settings must be chosen so that the POVM 
elements span the space of states, allowing reconstruction of any state in the space. 
In general there are optimal ways to choose both the settings and the amount of time 
spent measuring at each setting, which is sometimes referred to as optimal experiment 
design [7i [7i UHl IM] • 

Once the settings have been chosen, the measurements are performed, usually 
with a fixed time spent at each a. In this time, the number of events corresponding 
to each outcome 7 are recorded, giving a set of counts {n-yjc,. The probability that 
any state, represented in density matrix form as p, could produce the observed counts 
is then given by the likelihood function, 

C = p{{naj}\p), (13) 

where {ua-y} denotes the set of all measured counts at all settings. We discuss here 
the most common state reconstruction technique, known as maximum likelihood 
estimation, where the estimated p is determined by maximising the value of the 
likelihood function [76 1 177 ] . 

The exact form of C is dependent upon the system generating the state. With 
heralded state generation, the number of times that the state has been prepared for 
each setting, n„, is known so that for each herald signal a result is collected. The 
number of heralds limits the values of n-y so that C is most appropriately described 
by a multinomial distribution, 

^ = n^"(K}aK7' (14) 

where Cq, is a multinomial factor which accounts for the number of ways the observed 
n-y could have occured for each setting and Pa-f — tr{plla^). Since Ua is known for 
each setting, only N-y — 1 measured counts, or their corresponding probabilities, are 
independent variables in the maximum likelihood optimisation. This can be expressed 
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as the constraint equation ^ n^^ = Ua ■ This maximisation can be performed directly 
using an iterative technique [77] . 

If Ha is large and Paj is small, the multinomial form can be approximated by a 
product of Na{N^ — 1) Gaussian distributed variables, 



where aa^y = naPa^i^ ^Pay) ~ naPa-y for small pa-y- Typical current devices exhibit 
high loss, so that the vacuum component of any heralded state will contain a large 
population. While this might be expected to take us out of the regime of small Paj, 
we are still free to choose which outcome probability is the "dependent variable" for 
each measurement setting. Using the corresponding constraint equation, we therefore 
choose to eliminate the vacuum term, leaving only the lower probability terms which 
still satisfy the small p approximation. Expressing the likelihood function in this form 
ultimately allows one to express the problem in a weighted least-squares form which 
can be converted into a semi-definite program so that the well developed tools of 
convex optimisation can be employed [82] . 

A necessary part of this optimisation is knowledge of the POVM elements which 
model the particular outcomes observed. Since we wish to be able to reconstruct all 
photon number subspaces, we need measurements which access all photon number 
subspaces. In the case of a perfectly efficient, perfectly photon number resolving 
device, it can be precisely known which photon number subspace is being projected 
onto. In the presence of any loss, however, even photon number resolving detectors 
project onto multiple photon number subspaces since registering a single photon could 
be due to, for example, two photons where one is lost. The exact value of the loss 
in the measurement device specifies the mixture of photon number subspaces being 
projected onto and hence must be known if a state is to be accurately reconstructed. As 
an example, we will consider the device shown in figure[7]which provides tomographic 
measurements suitable for characterising the IIB(l) state. The detectors used are 
perfect APDs which have a POVM set of {|0)(0|,1 - |0)(0|} meaning that they 
either don't fire and therefore project onto the vacuum, or they fire in which case 
at least one photon was present but it is not known how many. This device is 
insensitive to coherences between components in different photon number subspaces. 
Fortunately, however, these coherences cannot affect the measurement statistics in our 
system, or indeed in any system comprised of imperfect state preparation, lossy linear 
optical interactions, and inefficient detection. This is true because all linear optical 
components conserve photon number by definition (i.e. they operate only within each 
photon number subspace), loss is always incoherent, and photon counting detectors 
project onto incoherent mixtures of components from within these different subspaces. 
A state reconstructed in this way therefore provides sufficient information to predict 
all measurement outcomes of any subsequent system acting on that state using linear 
optical interactions and photon counting measurement. 

For the system in figure [3 there are five possible outcomes: no clicks; a click at 
either of detector 1 (Dl) or detector 2 (D2); a click at detector 3 (D3); a click at Dl 
and D2; a click at Dl and D3 or D2 and D3. These are labelled n^^°"*^ where 7 = 1 
to 5. Losses in the system are modelled by beam splitters with refiectivities 771 to rjQ. 
Detector inefficiencies are incorporated into the losses directly before the detectors (775 
and tjq) [83| . The modes leading to detectors are labelled with mi to 7715 while the 
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loss modes are labelled l\o to Zeo for the "outer" modes and l\i to Zgi for the "inner" 
modes. Typically, the projectors for the detectors are back propagated and expressed 
in terms of projectors at the input of the measurement device and the r]i are set to 1. 
If losses are considered, however, one must back propagate both the click projectors 
for the APDs and identity projectors for each of the outer loss modes since they are 



not monitored. The output POVM set is given by the projectors 

n(°"*) = 1,^^ ^ ... ® 1,^^ |o)(oL3 ® |o)(o|,„4 ® loxoLs, (i6) 

n(°"*) = ® ... ® ® [(1 - |0)(0|„3) ® \mmA + lOXOLs ® (1 - \m\mA)\ 

«|0)(0L5, (17) 

ni°"*' = izio ® - ® i/eo ® |0)(0L3 ® |0)(0L4 ® (1 - |o)(o|„5) , (is) 
ni°"*^ = iho ® - ® ^ (1 - |0)(0L3) ^ (1 - |0)(0L4) ® |0)(0L5, (i9) 
n^°"*^ = i/io ® - ® i;e<, ® [(1 - |0)(0L3) ® + |o)(oL3 ^ (i - |o)(oL4)] 

^(l-lOXOLs). (20) 



By expanding the photon- number modes as |n)(n I = l/n!oT"|OXO|o", these ideal 
POVM elements can then be propagated backwards through the lossy circuit to 
determine the complete POVM elements for the overall measurement apparatus, 
which we represent by These back-propagated projectors, however, can be 

dramatically simplified using a priori knowledge of the system, that is, that the input 
loss modes In to l^i contain only the vacuum, and that the state p contains, for 
example, no more than two photons. The back-propagated projectors can therefore 
be conditioned using this information, which can be summarised by the following 
projector: 

P = (|0)(0| + |1X1| + |2)(2|)„i ® (|0)(0| + |1X1| + |2)(2|)„2 ® \m\iu ® - ^ |0)(0|j6,.(21) 
This allows us to calculate the final measurement POVM elements according to: 

= pn2;'p. (22) 

We note that when the transmissivities r/i to rj^ are pair-wise symmetric so that 

rji = ri2 and rj^ = 7^4, the situation is further simplified, because the total population 
in the loss modes is not a function of the phases 9, </> since every possible path 
to the detectors experiences the same loss regardless of the phase settings. This 
approximation corresponds to the standard bulk-optics scenario where loss does not 
typically depend on the occupied mode, for example, in a free-space polarization-based 
interferometer where loss is not normally polarization dependent. 

6. Conclusion 

Integrated photonics is one of several extremely promising technologies whose aims 
are to allow the exploration of the fundamental properties of quantum mechanics 
and the application of these phenomena to real tasks. Although much progress has 
been made in investigating the individual modules necessary to perform integrated 
experiments, there has been little work on managing the interfaces between them in 
order to achieve full integration. In particular, the effects of loss in quantum devices, 
which differ significantly from the effects of loss in classical devices, have been largely 
overlooked to date although they are now the single largest hurdle confronting the field. 
In this paper we have explored a modular approach to integration and enumerated the 
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requirements on each of state generation, manipulation, and detection in the context 
of quantum-enhanced metrology. Wc have argued that the control of modes is key, 
both in order to optimise each individual module, and to optimise the interfaces. 

We have discussed sources based on spontaneous four-wave mixing and the 
necessary considerations which enable pure, single photons to be generated with high 
heralding efficiencies in a transverse spatial mode that is well matched to existing 
photonic circuit technology. We have shown how managing the spectral mode of the 
generated photons is critical to maximising the heralding efficiency, demonstrating 
the trade off between this and state purity for a typically correlated source. We 
have also discussed how to construct programmable devices which are robust to 
fabrication imperfections and demonstrated a method that allows the characterisation 
of individual beam splitters, even when they are embedded within a complex circuit. 
Finally, we have enumerated some of the possibilities for integrated detectors and 
introduced a framework which describes how to perform state tomography of a 
heralded state, allowing reconstruction of all photon number sub-spaces rather than 
post-selecting on a particular photon number sub-space. This framework ultimately 
allows the utility of a state for quantum-enhanced metrology to be assessed through 
the quantum Cramer-Rao bound. 
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